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Investigations of the exactly solvable excitation spectra of 
two-electron quantum dots with a parabolic confinement, for 
different values of the parameter Rw expressing the relative 
magnitudes of the interelectron repulsion and the zero-point 
kinetic energy, reveal for large Rw a. ro- vibrational spectrum 
associated with a linear trimeric rigid molecule composed of 
the two electrons and the infinitely heavy confining dot. This 
spectrum transforms to that of a "floppy" molecule for smaller 
values of Rw- The conditional probability distribution calcu- 
lated for the exact two-electron wave functions allows for the 
identification of the ro- vibrational excitations as rotations and 
stretching/bending vibrations, and provides direct evidence 
pertaining to the formation of such molecules. 
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The behavior of three-body systems has been a con- 
tinuing subject of interest and a source of discoveries 
in various branches of physics, both in the classical and 
quantum regimes, with the moon-earth-sun system |l]] 
and helium-like atoms j^-Q] ( m the ground and excited 
states) being perhaps the best known examples. Further- 
more, insights gained through such investigations often 
provide the foundations for understanding the properties 
of systems with a larger number of interacting particles. 

Recently, analysis of the measured conductance || and 
differential capacitance ^ spectra of two-dimensional 
(2D) quantum dots (QD's,) created via voltage gates at 
semiconductor heterointerfaces, led to their naming (by 
analogy) as "artificial atoms". In particular, this anal- 
ogy refers to identification of regularities in the measure- 
ments which have been interpreted ^ along the lines of 
the electronic shell model (SM) of natural atoms, which 
is founded on the physical picture of electrons moving in 
a spherical central field including the averaged contribu- 
tion from electron-electron interactions. 

Motivated by the central role that spectroscopy played 
in the development of our undestanding of atomic struc- 
ture, we investigate in this paper the exactly solvable 
excitation spectrum of a two-electron (2e) parabolic QD 
as a prototypical three-body problem comprised of the 
two electrons (X's) and the (infinitely heavy) confining 
quantum dot (Y). Through probing of the structure 
of the exact wave functions with the use of the condi- 
tional probability distribution (CPD) ||], in conjunction 
with identification of regularities of the excitation spec- 
trum, we show that such a spectrum is characteristic of 
collective dynamics resulting from formation of a linear 
trimeric molecule XYX jl(|. In particular, we find that 
the excitation spectrum of the 2e QD exhibits for a weak 



parabolic confinement (i.e., small harmonic frequency 
luq) a well-developed, separable ro-vibrational pattern 
which is akin to the characteristic spectrum of natural 
"rigid" triatomic molecules (i.e., molecules with stretch- 
ing and bending vibrational frequencies higher than the 
rotational one). For stronger confinements (i.e., large 
ujq), the spectrum transforms to one characteristic of 
a "floppy" triatomic molecule, converging finally to the 
independent-particle picture associated with the circular 
central mean field of the QD. 

The Schrodinger equation for a 2e QD with a parabolic 
confinement of frequency loq, with the 2D Hamilto- 
nian given by H = Ei=i.2 P?/ 2m * + e 2 //-c|ri - r 2 | + 
0.5w* Uq J2i=i 2 r f ' wnere K an d rn* are, respectively, the 
dielectric constant and electron effective mass, is separa- 
ble in the center-of-mass (CM) and relative-motion (rm) 
coordinates Consequently, the energy eigenvalues 

may be written as EjsiM.nm = £™i +s" n (n, \m\), where 
= ^ (27V + \M\ + 1) with the N and M quantum 
numbers corresponding to the number of radial nodes in 
the CM wave function and M is the CM azimuthal quan- 
tum number, and e rm (n, |m|) are the eigenvalues of the 
one-dimensional Schrodinger equation 
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where Q(u)/^/u is the radial part of the rm wave func- 
tion Q(u)e lm0 /y/u with n being the number of radial 
nodes; u = |ui — u 2 | with = Ti/loV^ (i = 1,2) be- 
ing the electrons' coordinates in dimensionless units and 

= (h/m*LijQ) 1 ^ 2 , that is the spatial extent of the lowest- 
state wave function of a single electron. The so-called 
Wigner parameter Rw = (e 2 / kIq) / hujQ multiplying the 
Coulomb repulsion term expresses the relative strength 
of the Coulomb repulsion between two electrons sepa- 
rated by Iq and twice the zero-point kinetic energy of an 
electron moving in a harmonic confinement. 

Denoting the exact spatial wave function of the 2e QD 
by $jvM,nm(ui, u 2 ) (which is the product of the CM and 
rm wave functions) , and the spatial two-electron density 
by WjVM,nm(ui,u 2 ) = | $jvAf,nm(ui , u 2 ) | 2 , we define the 
usual pair-correlation function (PCF) as 



G(v) = 2tt J J <5(ui - u 2 - v)VF(ui, u 2 )duidu 2 , 

and the conditional probability distribution (CPD) for 
finding one electron at v given that the other is at Vo as, 

W(v, u 2 = v ) 
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FIG. 1. Spectra (bottom) and corresponding conditional probability distributions (CPD) and pair-correlation functions 
(PCF) shown respectively on the left and right of each of the subplots [labeled (a)-(d)], for a 2eQD with Rw = 200. For 
each excitation band, the quantum numbers (No, Mo, no,m) are given at the bottom with m — 0,1, 2, ... (the levels for m — 
and m = 1 are not resolved on the scale of the figure and appear as a thick line); only a few of the low lying rotational and 
vibrational states are shown, with the collective rovibrational behavior extending to higher excitations. The CPD's and PCF's 
are labeled with the quantum numbers of the corresponding levels. For the spectral rules governing the spectrum and for the 
definition and interpretation of the CPD's and PCF's, see the text. The solid dot in each of the CPD subplots denotes the point 
vo = (do,0), where do = 2.6 is half of the electron separation found in the PCF of the ground state (0,0,0,0). All distances 
x,y,v and do are in units of /ov2, and energies are in units of huJo/2. 
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where the M, N, n, m indices of W (and therefore of 
G and V) have been suppressed. Note that the exact 
electron densities are circularly symmetric. 

With the above, we solved for the 2e QD energy spec- 
tra and wave functions for values of Rw = 200, 20 and 
3. We discuss first the Rw = 200 case whose spec- 
trum and selected PCF's and CPD's are displayed in 
Fig. 1. As can be seen immediately, for such a large 
value of Rw, the spectrum of the 2e QD (bottom part 
of Fig. 1) exhibits the following three well-developed 
regularities: (I) for every band (No, Mo, no, m), with 
m = 0,1,2,..., while No,Mq and no are kept constant 
(in the following the subscript "zero" denotes a number 
that is held constant in a particular sequence), the en- 
ergy spacing between two adjacent levels m and m + 1 
increases linearly in proportion to 2m + 1; the bands 
(No, ±M , n , ±m) are degenerate. Note that the lev- 
els are spin singlet or triplet for m even or odd, respec- 
tively, (II) the bands (0,M o ,0,m) and (iV o ,0, 0, m) cor- 
respond to excitations of the center-of-mass motion with 
Mo and 2Ao vibrational quanta (phonons) of energy fwjo, 
respectively, and (III) the bottom levels of the bands 
(0,0, no, to) form a one-dimensional harmonic-oscillator 
spectrum (n + l/2)hw s . 

The above three "spectral rules" specify a well- 
developed and separable ro- vibrational spectrum exhibit- 
ing collective rotations, as well as stretching and bending 
vibrations |Q. Indeed, neglecting an overall constant 
term, the above rules can be summarized as, 

E NM ,nm = Cm 2 +(n+ l/2)huj s + (2N + \M\ + l)hu b , 

where the rotational constant C « 0.037, the phonon for 
the stretching vibration has an energy Tilu s ss 3.50, and 
the phonon for the bending vibration coincides with that 
of the CM motion, i.e., huj b = huj Q = 2 @|U (all en- 
ergies are given in dimensionless units of Tiujo/2). Note 
that the rotational energy is proportional to to 2 , as is 
appropriate for 2D rotations, unlike the case of natural 
triatomic molecules where the rotational energy has a 
term proportional to 1(1 + 1), I being the quantum num- 
ber associated with the 3D angular momentum. Observe 
also that the bending vibration can carry by itself an 
angular momentum TiM and thus the rotational angular 
momentum hm does not necessarily coincide with the 
total angular momentum h(M + to). 

Further insight into the collective character of the spec- 
trum displayed in Fig. 1 can be gained by examining 
the CPD's and PCF's associated with selected states of 
the rotational bands (No, Mo, no, m) (the CPD's are dis- 
played to the left of the PCF's; notice that the PCF's 
are always circularly symmetric). The band (0,0,0, to), 
being purely rotational with zero phonon excitations, can 
be designated as the "yrast" band, in analogy with the 
custo mary terminology from the spectroscopy of rotating 
nuclei (lj. 

In Fig. 1(a), we display the CPD's and PCF's for three 
specific states of the yrast band, i.e., the (0,0,0,0), the 



(0,0,0,3), and the (0,0,0,6). The corresponding PCF's 
are all alike and centered around 2do — 5.2, which im- 
plies that the two electrons keep apart from each other 
at a distance 2do- Due to the circular symmetry of the 
PCF's, however, one can only conclude that the two elec- 
trons are moving on a thin circular shell of radius do ■ To 
reveal the formation of an electron molecule, one needs to 
consider further the corresponding CPD's [plotted in the 
left column with v = (do, 0); the point v is denoted by 
a solid dot]. In fact, the CPD's demonstrate that the two 
electrons reside at all instances at diametrically opposite 
points, thus forming a linear molecule XYX with two 
equal bonds (X — Y and Y — X) of length do- In addition, 
one can see that all three CPD's are practically identical, 
in spite of the fact that the angular momentum changes 
from m = (lower subplot) to m = 6 (upper subplot). 
This behavior, namely the constancy of the bond lengths 
irrespective of the rotational energy, properly character- 
izes the electron molecule as a rigid rotor. 

Turning our attention away from the yrast band, we 
focus next on the bands (0, 0, 1, to) and (0, 0, 2, m), which 
are rotational bands built upon one- and two-phonon ex- 
citations of the stretching vibrational mode. We have 
verified that the PCF's and the CPD's corresponding 
to these bands share with the yrast band the property 
that they do not change (at least for the levels displayed 
in Fig. 1) as a function of m. Thus it is sufficient to 
study the bottom states, i.e., those with m = 0, (0,0,1,0) 
and (0,0,2,0), whose corresponding PCF's and CPD's are 
displayed in the lower and upper subplots of Fig. 1(b), 
respectively. The PCF's demonstrate the presence of in- 
ternal excitations with one and two nodes in the relative 
motion, but they yield no further information regarding 
the electron molecule. The CPD's, however, plotted here 
for v = (do,0) [the point v is kept the same for all 
subplots in Fig. 1] immediately reveal the presence of 
excitations (specified by the number of their nodes, i.e., 
here one or two) associated with the vibrational mode of 
the XYX molecule along the interelecton axis (namely, 
the stretching vibration). 

By examining the corresponding CPD's, one can fur- 
ther demonstrate that the two degenerate rotational 
bands (0, 2, 0, to) and (1, 0, 0, to) are built upon the lowest 
two-phonon excitations of the bending vibrational mode 
of the linear molecule XYX. Again, we have verified 
that it is sufficient to consider the two states at the bot- 
tom of the bands, namely the (1,0,0,0) [see lower sub- 
plot of Fig. 1(c)] and the (0,2,0,0) [see upper subplot 
of Fig. 1(c)]. It can be seen that both CPD's describe 
vibrational excitations of the XYX molecule which are 
perpendicular to the interelectron axis (namely, bending 
vibrations), with the one associated with the (1,0,0,0) 
level having one node and the one associated with the 
(0, 2, 0, 0) having no nodes (this is in agreement with the 
fact that the normal mode associated with the bending 
vibrations is related to the 2D harmonic-oscillator de- 
scribing the CM motion). We note that the correspond- 
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FIG. 2. Spectra (bottom) and CPD's (top) of the corre- 
sponding ground-state levels for 2eQD's with (a) Rw = 20, 
and (b) Rw = 3. At each of the bottom panels, the spectrum 
of the yrast band [i.e., (0,0, 0,m); m = 0,1,2,...] is shown 
on the left, and the lowest levels of the bands (0,0,no,m) 
for no = 0,1,2,..., forming a stretching vibrational spec- 
trum (i.e., with constant spacing) are displayed on the right. 
The solid dot in each of the CPD subplots denotes the point 
vo = (do, 0), where do is half of the electron separation found 
in the corresponding ground-state PCF's (not shown here). 
Energies in units of Tlujq/2 and distances in units of ZoV2; 
e.g., for GaAs material parameters (m* = 0.067m e , k = 12.9) 
and Rw = 3, one has HtJo = 1.2 meV and Iq\/2 = 43.54 nm. 



ing PCF's [see right column in Fig. 1(c)] fail to describe 
(in fact they are completely unrelated to) the bending 
vibrations; indeed they are identical to the ones associ- 
ated with the yrast band [Fig. 1(a)] which is devoid of 
any vibrational excitations. 

The CPD and PCF of the bottom level (i.e., with 
m = 0) of the rotational band (1,0, l,m), which is 
built upon more complicated phonon excitations of mixed 
bending and stretching character (not shown in Fig. 1), 
are displayed in Fig. 1(d). It is easily seen that the CPD 
represents a vibrational motion of the electron molecule 
both along the interelectron axis (one excited stretching- 
mode phonon) and perpendicularly to this axis (two ex- 
cited bending-mode phonons). In fact, the CPD in Fig. 
1(d) can be viewed as a composite made out of two CPD's 
shown previously, one in the lower subplot of Fig. 1(b) 
and the other in the lower subplot of Fig. 1(c). Return- 
ing to Fig. 1(d), one can see again that, in contrast to 
the CPD which enables detailed probing of the excitation 
spectrum, the information which may be extracted from 
the corresponding PCF is rather limited. 



The rigidity of the electron molecule, which is so well 
established for Rw = 200, will naturally weaken as the 
parameter Rw decreases and the XYX molecule will start 
exhibiting an increasing degree of "floppiness" . Such 
floppiness can be best observed in the yrast band, which, 
beginning with the higher levels, will gradually deviate 
from the spectral rule (I) discussed above, and eventually 
it will become unrecognizable as a rotational band. This 
is illustrated in the lower subplot of Fig. 2(a) which dis- 
plays the yrast band for Rw = 20. Specifically, one can 
see that only the lowest four levels honor approximately 
rule (I), the higher ones tending to develop a constant 
energy spacing between adjacent levels [this spacing con- 
verges slowly to the energy spacing Tlujq (i.e., to the value 
2 in dimensionless units) of the parabolic confinement]. 
In the case Rw = 3, one can hardly identify any rota- 
tional sequence in the levels of the yrast band [plotted 
at the bottom subplot of Fig. 2(b)]. Indeed, although 
the energy spacing between the second and the third lev- 
els is larger than that between the first and the second 
levels (but with a ratio substantially different than 3/1), 
the spacing between higher levels approaches quickly the 
value 2 of the external confinement. 

However, in spite of the floppiness exhibited by the ex- 
citation spectra in Fig. 2, the (singlet) ground-state of 
the 2e QD for both Rw = 20 and Rw = 3 drastically 
deviates from the Is 2 closed-shell orbital configuration 
expected from the independent-particle picture. Rather, 
as demonstrated by the corresponding CPD's [top sub- 
plots in Fig. 2], in both these cases of smaller Rw's, the 
ground state is still associated with formation of rather 
well-developed XYX electron molecules, but with pro- 
gressively smaller bond lengths. Finally, we remark that 
the stretching vibrations are more robust and tend to 
better preserve a constant spacing between the bottom 
levels of the bands (0, 0, no, to) [these levels were grouped 
in a vibrational band (0, 0, n, 0) and are plotted on the 
right-hand-side of the lower subplots in Fig. 2]. 

The remarkable emergence of ro-vibrational excita- 
tions for parabolically confined 2eQD's, under magnetic- 
field-free conditions, provides direct evidence for the for- 
mation of electron molecules in QD's, with their rigid- 
ity controlled by the parameter Rw- Such electron 
molecules and associated collective excitation spectra are 
general properties [fj5[-|l7f of QD's (with greater spectral 
complexity in many-electron QD's), whose observations 
(and manipulations through controlled pinning of the col- 
lective rotations [15(b)]) form outstanding experimental 
challenges. 
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